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Abstract 

First data on coherent threshold 7r electroproduction from the deuteron taken by the Al Collaboration at the Mainz Microtron 
MAMI are presented. At a four-momentum transfer of q 2 = —0.1 GeV 2 /c 2 the full solid angle was covered up to a center- 
of-mass energy of 4 MeV above threshold. By means of a Rosenbluth separation the longitudinal threshold s wave multipole 
and an upper limit for the transverse threshold s wave multipole could be extracted and compared to predictions of Heavy 
Baryon Chiral Perturbation Theory. 
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1. Introduction 

The electroproduction of neutral pions is one of the 
significant testing grounds of Chiral Perturbation The- 
ory (ChPTh, see [1,2] for an overview of the field). 
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The data of the production from the proton [3,4] are 
consistent with calculations of Heavy Baryon Chiral 
Perturbation Theory [5] and are in reasonable agree- 
ment with this approach for the free nucleon case. For 
this comparison still six low energy constants, intro- 
duced into ChPTh by renormalizing counter terms, 
had to be adjusted to the data. However, for the free 
neutron a strong prediction without further freedom 
can be deduced from the proton production amplitude. 
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Unfortunately, no free neutron target exists and one 
has to cope as normal with the problem of a target 
bound in a complex nucleus. Therefore, a model is 
needed to separate the elementary production ampli- 
tudes and the nuclear effects introducing theoretical 
uncertainties. The most evident choice for investi- 
gating the neutron bound in a nucleus would be the 
quasi-free reaction. The momenta of the recoiling nu- 
cleons and the initial Fermi momenta are, however, 
of the same order and, therefore, the dynamics of the 
full reaction are not well under control. On the other 
hand, assuming that the proton production amplitude 
is known and that a reasonable model for the sim- 
plest nucleus is at hand today, the coherent production 
from the deuteron can be used. In this case, the initial 
and final Fermi motion can be included in the model 
description of the deuteron via structure form factors. 

For the photoproduction case, a first coherent 
threshold measurement of d(j, w°)d at SAL [6] ex- 
tracted the transverse threshold s wave multipole 
| I in good agreement with the predictions of a 
fourth order calculation in the framework of ChPTh 
[7]. In photoproduction the reaction is identified by 
the detection of the decay photons of the neutral 
pion, and the small incoherent break up channel has 
to be subtracted by model dependent assumptions. 
By contrast in an electroproduction experiment the 
recoil deuteron has to be detected due to large back- 
ground of hard and soft photons, thus there is no 
incoherent contribution. The low kinetic energy of 
the heavy recoil deuteron (25 MeV < T d < 35 MeV) 
is the biggest challenge of such an experiment, the 
detection efficiencies and multiple scattering effects 
have to be measured at each state of the experiment. 

On the other hand, the detection of the deuteron 
with a high resolution spectrometer allows for a good 
invariant mass resolution close to threshold. The fo- 
cusing due to the Lorentz boost permits to cover the 
full solid angle in the center-of-mass system up to 
4 MeV above threshold if one uses the large solid an- 
gle magnetic spectrometer A of of the A 1 Collabora- 
tion with Afi = 21 msr. 

2. Kinematics 

In Born approximation the virtual photon is defined 
by the four-vectors of incident electron e = (E, k) 
and scattered electron e' = (E' , k') as q = e — e' = 




Fig. 1. Definition of angles. 

(io,q) (see figure 1). Denoting the variables in the 
photon-deuteron center-of-mass frame with an aster- 
isk the cross section can be written as 



da 



da 



dE'dQ'dfl* dfi* 
with the virtual photon flux 



(1) 



2tt 2 E -q 2 1 - e 



(2) 



and the equivalent photon energy fc 7 = (W 2 — 
m 2 d )/2md- The transverse and longitudinal degree of 
polarization of the virtual photon is given by 



f = ,l- 2 ^tan^ 
q z 2 



(3) 
(4) 



For the coherent production of pseudo scalar mesons 
from a spin 1 target the unpolarized differential cross 
section can be separated similarly to the production 
off a proton into four structure functions (see e.g. 
[8,9]): 



dnt 



+ v /2e L (l-e)-/ T L(^)-cosC 

+ e-/ T T(^)-cos2C. (5) 

In principle, the extraction of the transverse- 
longitudinal /tl and transverse-transverse interfer- 
ence structure functions ,/tt is possible by a mea- 
surement of the </>* dependence of the differential 
cross section. In the presented experiment, how- 
ever, the angular resolution of the detection of the 



2 



Table 1 

Kinematical settings. The four-momentum transfer for all 
settings is q 2 = -0.1 GeV 2 /c 2 . 



e 


Eo 
[MeV] 


E' 
[MeV] 


Pd 
[MeV/c] 


Be 


e d 


0.854 


720 


554 


339 


29.00° 


48.69° 


0.590 


435 


269 


339 


55.10° 


38.06° 


0.364 


345 


178 


339 


79.22° 


29.40° 



low energy deuterons was not sufficient to allow an 
extraction of these small structure functions. 

The choice of the four-momentum transfer of q 2 = 
—0.1 GcV 2 /c 2 as a lower limit was dictated by the 
detection efficiency of the deuteron as will be dis- 
cussed in the next section. For this four-momentum 
transfer we chose three different values for the virtual 
photon polarization e with the largest possible spread 
for a Rosenbluth separation. Table 1 summarizes the 
kinematical settings. 

3. Experimental Setup and Analysis 

3.1. Particle Detection and Efficiencies 

The experiment was performed at the three spec- 
trometer setup of the Al Collaboration at the MAMI 
accelerator (see ref. [10] for a detailed description 
of the setup). For the electron detection spectrometer 
B, a clamshell dipole spectrometer with an angular 
acceptance of 5.6 msr at a momentum resolution of 
Ap/p = 10~ 4 was used. For the deuteron detection 
spectrometer A with a large solid angle of 21 msr was 
chosen. 

A high power liquid deuteron target was used at 
luminosities of 15 MHz//xb, limited by the current in 
the drift chambers of the deuteron spectrometer. Spe- 
cial care had to be taken to minimize the pathlength 
of the deuterons in the target material. A long narrow 
cell of 4.8 cm length and 1 cm width was shifted to- 
wards the side of the scattered electrons. In this way 
an average pathlength of 3 mm for the deuterons in 
the target material could be achieved. In addition the 
beam had to be moved by a fast magnet across the 
target area in a time scale of several kHz to avoid 
boiling of the liquid deuterium. The wall of the target 
was built of a 10 fxm Havar foil. 

To further reduce multiple scattering and energy 
loss the vacuum of the scattering chamber was con- 



nected with the vacuum of the deuteron spectrome- 
ter. For the electron detection a focal plane detector 
consisting of 4 layers of vertical drift chambers for 
spatial resolution and two layers of scintillators for 
coincidence trigger and time of flight measurement 
was used. A halocarbon gas Cerenkov detector with 
an efficency of 99.8% was used for the separation of 
electrons and charged pions dominantly produced off 
the target walls. For the deuteron detection, only one 
layer of scintillators after four layers of vertical drift 
chambers could be used, since the deuterons were al- 
ready stopped after a pathlength of a few millimeters 
in the first scintillator layer. 

The large energy loss and the deuteron loss due to 
nuclear reactions of the low energy deuterons made 
it necessary to monitor the efficiency of the deuteron 
detection very carefully. In order to estimate these ef- 
fects we used the standard formulas for the energy 
loss [11] and the calculations of [12] for the contribu- 
tions of nuclear reactions. Since these calculations are 
only valid with large restrictions in our energy range, 
we checked and calibrated in addition the deuteron 
detection efficiency by a coincidence measurement of 
the scattered electron and the deuteron in the elastic 
d(e,e'd) reaction. We used several settings of the elas- 
tic line at different positions of the focal plane and 
compared the results with the known cross section in 
the parameterization of [ 1 3] . 

As expected from extensive simulations on the 
computer, we had to apply correction factors of the 
order of up to 1.5 for the lowest deuteron energies. 
Figure 2 shows this correction for the elastic mea- 
surements in comparison with the calculated elastic 
cross section. 

3.2. Reaction Identification 

The reaction d(e, e'd)ir was identified in two 
steps. First, the coincidence time was determined by 
measuring the time of flight of the deuteron and the 
electron and correcting it for the reconstructed path 
length inside the spectrometers and the measured 
momenta. The coincidence time resolution of 3.3 ns 
FWHM was limited by the uncertainty in the flight 
path reconstruction because of the large multiple 
scattering of the low energy deuterons. 

After the timing cut, the missing mass was cal- 
culated from the measured four-vectors of incident 
electron e, scattered electron e', and initial and final 
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Fig. 2. The elastic cross section of the d(e,e'd) reaction 
at an incident energy of Eo = 420 MeV. The solid line 
shows the known cross section calculated with the param- 
eterization of [13]. The open circles show the coincidence 
measurement without efficiency correction, the solid cir- 
cles show the measured cross section after full efficiency 
correction. 

deuteron d and d' by w^iss = (e + d — e' — d') 2 . 
Figure 3 shows the distribution of the missing mass 
before and after background subtraction. As can be 
seen, the tt° can be identified with a missing mass 
resolution of 2.28MeV/c 2 FWHM, which is again 
determined by the multiple scattering of the low en- 
ergy deuterons in the wire chambers. 

The effect of the radiative tail at higher missing 
masses was corrected by a simulation using the stan- 
dard formulas of [14] in the peaking approximation. 

Since the recoil deuteron is measured directly, no 
model dependent correction for the deuteron break up 
is necessary. 

3.3. Systematic errors 

The error of this experiment is dominated by the 
systematic errors. As known for threshold measure- 
ments the calibration of the measured electron mo- 
mentum which is almost proportional to the center- 
of-mass energy causes the largest contribution. 
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Fig. 3. The missing mass distribution with a cut on the 
coincidence time (total) and with a cut on a timing win- 
dow shifted by 5 ns besides the coincidence time peak for 
the background (dark gray). The lower picture shows the 
background subtracted missing mass spectrum. The radia- 
tive tail appears at higher masses. 

Table 2 

Statistical and systematical errors of the total cross section 



Setting 


e = 


0.854 


e = 


0.590 


e — 


0.364 


AW 


stat. 


sys. 


stat. 


sys. 


stat. 


sys. 


[MeV] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


0.5 


7.5 


22.7 


6.9 


20.3 


7.1 


22.9 


1.5 


3.7 


12.7 


3.2 


9.1 


3.8 


10.4 


2.5 


2.9 


6.0 


2.8 


3.3 


3.0 


4.1 


3.5 


2.8 


3.6 


3.0 


3.7 


3.2 


4.8 



The systematic error caused by the efficiency cor- 
rection of the deuteron detection could be checked by 
our elastic calibration measurements and was deter- 
mined to be 1.7%. Further the detector efficiencies, 
contributions of cuts and phase space integration, er- 
rors of the luminosity summation, and condensation 
on the cold target walls were taken into account. 

The total systematical errors, together with the sta- 
tistical errors, are compiled in table 2. 
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Fig. 4. The differential cross section for e — 0.590. 

4. Discussion of the Results 

The measured differential cross section is compiled 
in table 3. For one setting (e = 0.590), figure 4 shows 
the data points, including the combined statistical and 
systematical errors. As expected, only angular struc- 
tures up to O(cos 2 9* ) appear and justify a fit with 
the assumption of only s and p waves contributing to 
the cross section. Figure 5 shows the total cross sec- 
tion and the result of a \ 2 fit to the data with the as- 
sumption of a constant s wave amplitude and p wave 
amplitudes rising linearly^] with p*. At the precision 
of these data, no cusp effects of the opening deuteron 
break up threshold at 2.2 MeV above threshold can 
be observed. 

At present, only threshold calculations are available 
in the framework of ChPTh [15] which provide a 
prediction for the s wave multipoles. From the fit we 
can extract the s wave multipoles at threshold through 
the reduced threshold s wave cross section 



_ 3 k* da 



\E d \ 2 + e L \L d \ 2 



(6) 



3 In this energy range the difference between a lineal' behavior 
as found empirically [6] and the theoretically predicted behavior 
~ p* /i±f can be neglected. 
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Fig. 5. The total cross section for three different values of 
the photon polarization e. The lines show the result of a 
least squares fit with the assumption of only p and s waves 
contributing to the cross section near threshold. 

Since 6l ~ 9e the longitudinal amplitude L d con- 
tributes with a much higher weight to the s wave 
cross section than E d and only an upper limit of one 
standard deviation can be extracted for From this 
analysis one gets 

\E d \ <0.42- IQ-^/m^ 

\L d \ = (0.50 ±0.11) • 10~ 3 /m v 

The classical Rosenbluth method, i.e. the total cross 
section plotted against the photon polarization e, rep- 
resents a consistency check (Figure 6). For each en- 
ergy bin a straight line fit was performed to separate 
/t by its offset and /l by its slope. At threshold, 
only |Ld| 2 contributes to the longitudinal cross sec- 
tion and can be determined by an extrapolation of 
/l to the threshold point. By this technique \L d \ = 
(0.47 ± 0.18) ■ lO^/m^ is extracted in good agree- 
ment with the previous analysis. Again, the kinemati- 
cally suppressed multipole E d cannot be determined. 

Figure 7 shows the extracted s wave multipoles in 
comparison with the prediction of Chiral Perturba- 
tion Theory [15]. These calculations are performed 
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Fig. 6. Rosenbluth plot with straight line fits including one 
standard deviation error band. 



to third order in the chiral expansion and are shifted 
to reproduce the result of the fourth order calcula- 
tion at the photon point [7]. The solid line shows the 
full calculation, for the dashed lines the calculated 
free neutron amplitude was varied by ±10~ 3 /m 7I to 
indicate the sensitivity of the calculation to this am- 
plitude. The dash-dotted line shows the calculation 
without two body currents, i.e. only the amplitudes of 
the free nucleons folded with the pertinent deuteron 
form factors are included and no pion exchange be- 
tween the two nucleons was taken into account. For 
this picture, we assumed the sign for id to be the 
same as calculated in ChPTh. 

As stated in [6], their transverse multipole is 
already at the photon point AE d = 0.35 • l0~ 3 /m n 
above the ChPTh prediction [7]. A shift by this 
amount would make our result for consistent with 
the calculations of [15]. 



+ Ewald et al. 
4 Bergstrom et al. 



- single Scattering 




0.05 
-q 2 [GeV 2 /c 2 ] 



0.1 




0.05 
-q 2 [GeV 2 /c 2 ] 



0.1 



Fig. 7. The extracted s wave multipoles (circles) in com- 
parison with the prediction of Chiral Perturbation Theory 
[15]. The photon point of [6] is plotted as a square. For 
the explanation of the curves see text. 



5. Summary 

A first measurement of the coherent threshold elec- 
troproduction of neutral pions off the deuteron was 
performed and the differential cross section could be 
determined up to 4MeV above threshold for three 
different values of the photon polarization. The abso- 
lute values of the s wave multipoles were extracted 
and compared to the predictions of ChPTh. Although 
smaller than expected, our extracted upper limit for 
\E A \ is consistent with ChPTh calculations. |i d | is 
overestimated by the theory by a factor of 2. More 
dramatic is this discrepancy in terms of cross sec- 
tions: The reduced threshold s wave cross section an 
is one order of magnitude smaller than expected. 

Both, in photo- and in electroproduction the mea- 
sured differential cross sections allow tests of predic- 
tions for the contributing p waves which are not yet 
calculated in the framework of ChPTh. 
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Table 3 

The differential cross section. The cross section is inte- 
grated over the complete accepted out-of-plane angle 0* . 

da/dSl* [nb/sr] 
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-0 


,083 


55 


.06 ± 


7.23 


22.19 ±4.95 


22. 


,52=b4. 


.40 


0. 


,083 


41 


.72 ± 


6.70 


31.30 ±5.77 


32. 


23 ±5. 


.26 





,250 


41 


.41 ± 


7.06 


45.12 ±6.52 


25 


.SI ±4. 


.89 





,417 


26 


.78 ± 


6.49 


30.65 ±5.92 


23. 


35 ±4, 


.76 





,583 


24 


.05 ± 


6.84 


32.06 ±6.23 


27. 


.31 zb4. 


.74 


0. 


.750 


29 


.06 ± 


8.52 


29.70 ± 6.98 


28. 


.Ilzb4, 


.99 





,917 


16 


.78 ±10.93 


23.53 ±8.55 


23. 


.70 ±5. 


.06 
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